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Abstract 

I review recent theoretical progress in the resummation of small x con- 
tributions to the evolution of parton distributions, in view of its poten- 
tial significance for accurate phenomenology at future colliders. I show 
that a consistent perturbative resummation of collinear and energy logs 
is now possible, and necessary if one wishes to use recent NNLO results 
in the HERA kinematic region. 
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I review recent theoretical progress in the resummation of small x contributions to 
the evolution of parton distributions, in view of its potential significance for accurate 
phenomenology at future colliders. I show that a consistent perturbative resummation 
of collinear and energy logs is now possible, and necessary if one wishes to use recent 
NNLO results in the HERA kinematic region. 

1. The dangerous success of NLO calculations 

There is a general attitude in the collider physics community that small x resum- 
mation is impossible to understand, but there is no need to worry since the data can 
be described ver y wel l without it. It is certainly true that (for instance) standard 
global parton fits UJ^I are based on fixed next-to-leading (NLO) order calculations, 
and that they manage to describe the data very accurately, in particular through- 
out the HERA kinematic region. The lack of experimental evidence for higher-order 
perturbative corrections at small x, despite theoretical arguments which suggest 
that they should be large, has been the mot ivation for a substantial amount of 
theoretical effort over the last decade. This activity has become less of an 

academic exerc ise since the recent determination of the full set of NNLO splitting 
functions^^^and of the full set of 0(a s ) perturbative corrections to deep-in elastic 
scattering. Indeed, parton fits that include these NNLO terms appear 

to be unstable at small x: when x 5s 10 -3 the difference in results e.g. for the 
P gg splitting function or for the gluon-dominated structure function Fl when go- 
ing from NNLO and NLO starts becoming as large or larger than the difference 
between NLO and LO. Interestingly, it also turns out that at small x there is a sig- 
nificant difference between the full NNLO result and the nominally leading small x 
contributions: subleading terms are crucial. This means that, whatever the reason 
for the success of NLO fits, small x resummation is mandatory beyond NLLO, at 
least for x < I0~ 3 . 

Thanks to recent theoretical progress, there is now a consistent theory which 
enables small x resummation, and which is close to being amenable to realistic phe- 
nomenology. This theory requires several ingredients, which were developed by vari- 
ous people over the last decade, and have recently led to two approaches (ABfE1"E3 
and CCSS ^B-IMI) which incorporate similar basic principles, and which arrive at 





1 



February 2, 2008 4:8 WSPC/Guidclines-IJMPA beijing 





1 2 3 4 5 



1 2 3 4 5 



Fig. 1. Dual expansions of the BFKL kernel \ and the GLAP anomalous dimension 7. 

stable and consistent resummed results. Here we will mainly review the ABF ap- 
proach, while also comparing with the results of the CCSS approach and briefly 
discussing differences between the two approaches. 

The main problems in small x resummation of the standard GLAP evolution 
equations is first, that leading logarithmic resummation corrections are to be too 
large — they would lead to a powerlike small x rise of splitting functions and thus 
of structure functions which is incompatible with the data, 11 1 and second, that 
they are perturbatively unstable — the next-to-leading logarithmic resummation 
corrections are larger than the leading ones. 

These problems are solved thanks to three main ingredients: duality l ^4 | 4 | 5 | 
which makes the joint resummation of collinear and small x logs possible, running 
coupling small x resummation^Eland factorization ^, which s oftens the resummed 
small x behaviour, and gluon exchange symmetry | H | 14 | 10j which stabilizes the 
resummed perturbative expansion. 



2. Duality 

As is well known, singlet splitting functions at small x receive contributions of the 
form a s ±(a s In ~) n (c° + a sCpj + • ■ ■) to all perturbative orders. These contributions 
can be extracted from the BFKL equation, which for a parton distribution G(x, Q 2 ) 
(in fact, an eigenvector of the singlet evolution matrix), and its Mellin transforms 



G(N,t) 



d£e- N Z G&t) 



G(£, M) = I dte- m G{£,t) 



takes the form 



±G($, M) = X (M, a.) G(t M); £=-, 
at; X 



(1) 

(2) 
(3) 
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Fig. 2. Various expansions of the dual \ kernel, determined with n-f = 4, a s = 0.2. Note that 
the associate dual Eq. JSJ y(N) anomalous dimension is simply the inverse function of x(-^0i as 
per the axis labelling. 



of a standard renormalization-group (or GLAP) equation, 

±G(N,t)= 1 (N,a s )G(N,t); t=% (4) 
at n z 

but with the roles of the variables x and Q 2 interchanged. Upon Mellin transfer- 
mation, i^" 1 «-> (ln£)* and ^t^ 1 «- (ln£)*. 

Whereas the way to extract the coefficients of the leading singular (%■) contri- 
butions to the anomalous dimension ^(N) Eq. JU from the (BFKL) kernel x(M) 
of Eq. ll jfll has been known since a long timej^ it has been realized only more 
recently laililfll that in fact, up to higher twist corrections, the BFKL and GLAP 
equations admit the same solution, provided the boundary conditions are suitably 
matched, and the corresponding kernels satisfy the duality relations 

xh(N,a s ),a s ) =N (5) 
1 ( X (M,a s ),a s )=M. (6) 

These relations hold at fixed coupling, whereas in the running case they are cor- 
rected by terms which may be computed order by order in perturbation theory. 

Duality maps the expansion of 7 in powers of a s at fixed N (in Fig. 1 leading 
70, next-to-leading 71 etc.) into the expansion of \ in powers of of a s at fixed a s /M 
(in Fig. 1 leading Xs, next-to-leading Xss etc.): hence 70 or Xs sum leading logs of 
Q 2 (collinear logs) while 7 S or xo sum leading logs of ~ (energy logs) and so on. 
A joint (double-leading) resummation can be constructed by simply combining the 
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Fig. 3. Right: the singular running coupling correction to the NLO DL \ and it resummation. 
Left: running coupling resummation of the anomalous dimension 7 associated to a quadratic \ 
kernel. 



two expansions and subtracting the terms which are in common: so the LO DL 
anomalous dimension is 7dl,lo(A) = Oi s jo(N) + 7 S (a s /N) — and so on. A 
nontrivial property of the DL expansion is that the dual of LO DL \ is LO DL 7 
(up to subleading terms), and similarly at next-to-leading order and so on. 

The DL expansion allows one to understand and cure the notorious problem of 
the large size of subleading xi i n comparison to leading ^O: seen in Fig. 2, where 
BFKL LO denotes a s xo and BFKL NLO denotes a s Xo + c? s X- Indeed, momentum 
conservation implies that to all perturbative orders 7(1) = 0. By duality this also 
implies x(0) = 1. But x is a polynomial in a s , so x(0) = 1 means that as M — > 
X must behave as 

/ „ r s a a a 2 a 3 

up to subleading corrections. Hence the expansion of x in powers of a s has alternating- 
sign poles at M = 0, which are resummed into Xs and thus absent in double- leading 
X- Indeed, Fig. 2 also displays the DL LO and NLO curves, which agree respec- 
tively with LO and NLO GLAP when M £ 0, and with LO and NLO BFKL when 
M ~ h. Hence, the DL expansion is stable for all M k. 



3. Running coupling 

The running of the coupling a{t) — a M [l — /3oa^i+- . .] is a leading log Q 2 , but next- 
to-leading log i effect. As a consequence, beyond LLa; the fixed-coupling duality 
relations (JSJ get corrected: for instance "f ss determined from xo an d Xi according to 
Eq. O must be supplemented by a running-coupling correction A7 SS = — (3o ^°^\2 ■ 



1 



February 2, 2008 4:8 WSPC/Guidclines-IJMPA beijing 




Fig. 4. Symmetrization of the \ kernel of Fig. 2. 

One can view these running coupling corrections as a contribution to an "effective" 
X which then respects the duality relation JS} even at the running coupling level. 

The good thing about these running coupling correction is that one can show 
that through their inclusion duality between the BFKL and GLAP equations holds 
to all perturbative orders. This means that the running-coupling BFKL equation, 
just like the GLAP equation, admits a factorized solution, whose scale dependence 
is determined by the kernel independent of the boundary condition. El The bad thing 
is that the associate effective x is singular at the minimum of the fixed-coupling 
LO %, as shown in Fig. 3. This singularity implies that the associate splitting 
function grows as a power of £ in comparison to the leading-order one as x — ► 0: at 
N k LO, AP sk {a s ,x)/P(a s ,x) ~ {(3 a s k - These singularities must therefore be 
resummed to all orders if one wishes the x — > limit to be stable. 

The resummation can be accomplished in an asymptotic expansion, it leads 
to the nonsingular result of Fig. 3, and it has a further interesting consequence. 
Namely, if x nas a minimum, duality implies that 7 has a cut: e.g. if x is quadratic, 
then 7 has a square-root cut. After running coupling resummation, however, the 
cut in 7 is replaced by a simple pole (Fig. 3, right). Interestingly, the pole is located 
on the real axis to the left of the cut, as seen in Fig. 3 where a cut at TV ~ 0.5 
leads to a pole at N ~ 0.2. Because the location No of the rightmost singularity 
in 7 implies a small x behaviour of parton distributions ~ x~ N ° this means that 
running coupling corrections considerably soften the resummed small x behaviour, 
a result first obtained in Ref. El 



4. Gluon exchange symmetry 

The double- leading perturbative expansion of x(M) is stable for M ^ 5, but still 
unstable in the vicinity of M = 1, where it has alternating-sign poles. In particular, 
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Fig. 5. Exact LO and NLO and approximate NNLO BFKL kernels with a s = 0.2. Note the 
overall prefactor of a s is omitted (unlike in previous Figs. 2—4). 

it will lack a minimum at even perturbative orders. This is problematic for the 
running coupling resummation discussed in the previous section, which relies on 
the existence of a minimum. The instability at M = 1 can be understood •= U| and 
cured E on the basis of the observation that in fact, due to the symmetry of the 
underlying Feynman diagrams upon interchange of incoming and outgoing gluons, 
the BFKL kernel is symmetric about M = i, i.e. x{M) = x(l~^f). This symmetry 
is broken by the DIS choice of kinematical variables, which treats asymmetrically 
the initial and final parton virtualities /i 2 and Q 2 , but it can be restored by choosing 
e.g. x = where s is the center-of-mass energy of the partonic process. It is also 
broken by the asymmetric choice of scale in running of the coupling Q^( C? 2 ). 

These symmetry breaking effects are computable and can be undone: once 
the symmetry is restored, one can symmetrize the DL expansion of Xi thereby 
obtaining a kernel which is perturbatively stable and free of poles at both M = and 
M = 1. In fact, the momentum conservation constraints implies that symmetrized 
X is an entire function of M, and has a minimum to all orders. The symmetrized 
DL expansion of x which ensues is displayed in Fig. 4, where one sees that the 
powerful combination of duality and gluon symmetry leads to a stable expansion: 
the LO and NLO approximations are quite close. One can then revert to DIS 
variables (Fig. 4, right), thus obtaining the kernel which lead through duality and 
running coupling resummation to an anomalous dimension and splitting function. 
Interestingly, symmetrization implies a further softening of the kernel: the minimum 
of Xi hence the branch cut of the corresponding 7, is moved (for a s = 0.2) from 
No ~ 0.5 to Nq ~ 0.3. Note that the NLO resummed behaviour is now harder than 
the LO one. 

The combination of symmetrization and running-coupling duality can be ex- 
ploited to obtain powerful analytic results. For instance, one can use the knowledge 
of 7 up to NNLoP ^ I ^^ I to determine all singular contributions to NNLO Xj as 
shown in Fig. 6. This requires a treatment of running-coupling corrections up to 
NLO, and of various interference terms, which is fe asible if running coupling duality 
equations are solved in an operator approach. EH 

5. Results 

Various approximations to the splitting function are displayed in Fig. 6, determined 
with rif = in order to avoid problems related to the diagonalization of the splitting 
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Fig. 6. The GLAP splitting function compared to various resummed approximations. 

function matrix. The GLAP LO and NLO results are very close to each other and 
coincide at small x (their small difference is proportional to n/), but the NNLO 

3 

GLAP result is seen to be unstable at small x because of an unresummed ^f- In — 
term (with negative coefficient). The simple (LO) fixed-coupling double leading 
result ('DL fix' in Fig. 6) displays a dramatic rise ~ a; -0 ' 5 which is certainly not 
seen in the data. This is considerably softened by the running coupling resummation 
('DL run' in Fig. 6). A yet softer behaviour is obtained if the running coupling 
resummation is applied to the symmetrized result ('res LO' and 'res NLO' in Fig. 6), 
which also leads to a stable perturbative expansion. 

For compariso n the N LO curve (i.e. including NLO resummation of £ and t) 
obtained by CCSS E3E21 j s a i so shown. The CCSS approach starts from the BFKL 
equation, with the aim of determining the off-shell gluon density, not just the 
anomalous dimensions of parton distributions. This is then improved by enforc- 
ing symmetry, matching to GLAP (analogous to duality), and running coupling. 
The main advantages of this approach are its somewhat wider scope, and the fact 
that, being based on the BFKL equation, it allows for an exact treatment of the 
running coupling at small x (unlike the ABF approach, where it is determined 
through an asymptotic expansion). The main disadvantages are that, due to the 
latter feature, anomalous dimensions can only be obtained through a numerical de- 
convolution procedure, and due to the former feature the perturbative expansion of 
the anomalous dimension is not naturally organized as a leading log series (unlike 
the DL approach), which implies for instance that the full NLO nf = 4 result in 
this approach is not yet available. 

The NLO ABF and CCSS results turn out to be quite close: their difference can 
be taken as an estimate of the intrinsic uncertainty in the resummation procedure. 
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These results are now essentially ready for a phenomenological implementation, 
though this will require the further development of suitable tools, such as matched 
resummed coefficient functions (for which all the required theoretical information 
is available), and interpolation of the resummed results for their efficient numerical 
implementation. 
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